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Abstract. We investigate the propagation of gravitational waves in the context of fourth 
order gravity nonminimally coupled to a massive scalar field. Using the damping of the orbital 
period of coalescing stellar binary systems, we impose constraints on the free parameters of 
extended gravity models. In particular, we find that the variation of the orbital period is 
a function of three mass scales which depend on the free parameters of the model under 
consideration; we can constrain these mass scales from current observational data. 
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1 Introduction 

Our Universe appears spatially flat and undergoing a period of accelerated expansion. Several 
observational data probe this picture [1-6] but two unrevealed ingredients are needed in order 
to achieve this dynamical scenario, namely dark matter at galactic and extragalactic scales 
and dark energy at cosmological scales. The dynamical evolution of self-gravitating structures 
can be explained within Newtonian gravity, but a dark matter component is required in order 
to obtain agreement with observations [7]. 

Lately, models of extended gravity [8-10] have been considered as a viable theoreti¬ 
cal mechanism to explain cosmic acceleration and galactic rotation curves. In such models 
one extends only the geometric sector, without introducing any exotic matter. Such models 
may result from an effective theory of a quantum gravity formulation, which may contain 
additional contributions with respect to General Relativity, at galactic, extra-galactic and 
cosmological scales where, otherwise, large amounts of unknown dark components are re¬ 
quired. In the context of models of extended gravity , one may consider that gravitational 
interaction acts differently at different scales, whilst the robust results of General Relativity 
at local and solar system scales are preserved [11]. 

Models of fourth order gravity have been studied in the Newtonian limit (weak-field 
and small velocity), as well as in the Minkowskian limit [12], In the former one finds modi¬ 
fications of the gravitational potential, whilst in the latter one obtains massive gravitational 
wave modes [13]. The weak-field limit of such proposals have to be tested against realistic 
self-gravitating systems. Galactic rotation curves, stellar hydrodynamics and gravitational 
lensing appear natural candidates as test-bed experiments [14-19]. 

Corrections to the gravitational Lagrangian were already considered by several au¬ 
thors [20-26]). From a conceptual viewpoint, there is no reason a priori to restrict the 
gravitational Lagrangian to a linear function of the Ricci scalar minimally coupled to mat¬ 
ter [27]. In particular, one may consider the generalization of f(R ) models, where R is 
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the Ricci scalar, through generic functions containing curvature invariants such as the Ricci 
squared ( R a pR a P) or the Riemann squared {Rap-ydR 01 ^ 6 ), which however are not invariant 
due to the Gauss-Bonnet invariant R 2 — 4 R^R^ + R flu x a R^ l '^ a . Note that the same remark 
applies to the Weyl invariant Hence, one may add a (massive propagating) 

scalar field coupled to geometry; this leads to the scalar-tensor fourth order gravity. 

At this point, let us comment on the possibility that extended gravity theories can be 
plagued by pathologies, such as the appearance of ghosts (negative norm states), which could 
allow for negative possibilities and consequently violation of unitarity [28-33]. In particular, 
while standard General Relativity and the Gauss-Bonnet theory have the same field content, 
this is not the case for the f(R) gravity type and Weyl gravity. The former is safe, even 
though it does not improve the ultraviolet behaviour of the theory; f(R) gravity has just 
an extra scalar and can be ghost free in its Newtonian limit [34]. The latter has an extra 
pathological spin-2 field, which however causes no problem in the low-energy regime since 
the effects of higher order terms give rise to small corrections to General Relativity [34, 35] . 
In our analysis, we will consider an action motivated from noncommutative geometry within 
the class of a Weyl type gravity. Since, as we will discuss, this proposal will be considered in 
the spirit of an effective field theory, it is free from any pathologies. 

Our aim here is to obtain, in the framework of post-Minkowskian approximation (weak- 
field limit) of a scalar-tensor fourth order gravity , a general solution for the wave propagation 
of scalar and tensor modes. The analysis will be carried out in the context of stellar binary 
systems. More specifically, by exploiting recent astrophysical data on the variation of the 
orbital period of binary systems, we will constrain the free parameters, namely the three 
masses {my, thr, m that characterize the scales on which higher order terms generated by 
the models of extended gravity become relevant. 

The outline of this paper is the following. In Sec. 2, we give the action of a scalar-tensor 
fourth order gravity and write down the corresponding field equations, which we then solve 
in the presence of matter within the weak-field approximation. In Sec 3, we compute the 
gravitational wave emission from a quadrupole source, and then discuss the energy loss in 
Sec. 4. Using astrophysical results on the orbital period damping, we infer lower limits on 
the free parameters of scalar-tensor fourth order gravity models studied in the literature in 
Sec. 5 and Sec. 6. We summarize our conclusions in Sec. 7. 


2 Scalar-Tensor Fourth Order Gravity 


Consider the action 


5 



f(R,R a pR a P, 


4>) + uj(cj))(t>. a (j)’ a + xc } 


( 2 . 1 ) 


where / is an unspecified function of the Ricci scalar R, the curvature invariant R a pR a P = Y 
(where R^ u is the Ricci tensor), and the scalar field (f>. The £ m is the minimally coupled 
ordinary matter Lagrangian density, a; is a generic function of the scalar field (j), g is the 
determinant of the metric tensor and X = 87tG. We use the convention c = 1. 

In the metric approach, the field equations are obtained by varying the action (2.1) with 
respect to gg v , leading to 

IrR„u - / + 3 ^’% ^ - fRu- + 9^fR + 2 fyR^Rau - 2[fyR a ( ^ u)a 

+D[/y^] + [fyRagl’^g^ + = x T u" , ( 2 - 2 ) 
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where T fll/ = —^= S ^'^jJiv rn ' > is the the energy-momentum tensor of matter, f R = 

fy = and O = = - a ,a is the D’Alembertian operator. The convention for the Ricci 

tensor is Rp U = R a RrTU , while for the Riemann tensor we define R a = rjg + with the 
affinities being the Christoffel symbols of the metric: T(V = \g ,irJ {g a a,p + gpa,a ~ gap,a)- The 
adopted signature is (-1-). The trace of the field equations (2.2) reads 

f R R + 2 f Y R a pR a ? - 2/ + D[3 f R + f Y R] + 2 [f Y R a %p - w(0)0 ;Q 0;“ = XT , (2.3) 

where T = T a a is the trace of energy-momentum tensor. Varying the action (2.1) with 
respect to scalar field 0 we get 


2w(0)[H0 + u^(0)0; a 0 ;a - fd> = o , (2.4) 

where w*(0) = - and U = %■ 

We will analyze the field equations within the weak-field approximation in a Minkowski 
background rj fiu \ 

g^u ~ Vfiu + h^ v , 0 ~ 0 (o) + p . (2.5) 

We develop the function / as 

f(R, R a pR af) , 0) ~ fn( 0,0, 0 (o) ) R + /^(° 2 0 > (0) ) R 2 + M ° ; Q^ (0) ) ^ _ 0 (o) )2 

+f R 40 , 0 ,<j> W )R(<j> - cj)W) + fy(0,0,^)R a pR a P ; ( 2 . 6 ) 


any other possible contribution to / is negligible [36-38]. The field equations (2.2), (2.3) and 
(2.4) then read 


(□r; + m Y 2 )Rfu> 


m R 2 - my 2 2 
■imn‘ ^ 





m^r^rny^ \ 
6 m R 2 v ) 


- my 2 fR4>{ 0,0,0 (o) ) ( d^ - g^U^P = m Y 2 V , 
(□r; + rn R )R - 3 m R f R(j> { 0,0, 0 (o) ) \J v p = -m R V T , 

(□rH ~ m^p - f R <p(0,0, (j)^) R = 0 , 


R 


(2.7) 


where IZ0 is the D’Alambertian operator in flat space and have made the definitions 


2 = _ fd Q,Q,0 (o) ) 

mR 3/hr( 0, 0, 0(°)) + 2/y (0,0, 0(°)) ’ 

2 = f R ( O,O,0 (o) ) 

^ /y(0,0,0W) ’ 

2 f H (0, 0,0(°)) 

20,(0(01) ■ 


( 2 . 8 ) 


The geometric quantities R pv and R are evaluated to the first order with respect to the 
perturbation h^. Note that for simplicity 1 we set f R ( 0,0, 0(°)) = 1 and u (0 (o) ) = 1/2. 
The Ricci tensor in Eq. (2.7), in the weak-field limit, reads 


R^v = h c 


in, v)o 


~ f-l?y kl[l,v ~ h 


,fu/ > 


(2.9) 


1 We can define a new gravitational constant: X 

/^(O,O,0°)/fl(O,O,0(°)). 


X f R (0,0,<t> w ) and f R4 ,{ 0,0,/ 
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where h = h a a . Using the harmonic gauge condition g pa T a p(J = 0 we have h Ato -’ <J —1/2 h )At = 
0, hence the Ricci tensor becomes Equation (2.7) then reads 


(I——1?7 T TTlY 


m /; 2 - my 2 q2 


L 3m^^ + V 

+ 2 my 2 /^(O, 0, </> (0) ) (<9 2 v - 

"> ^ _ _ o . /n\ . _ 


/ my 2 | mR 2 + 2 my 2 rn N 

v ^“ + 6 ^ ^ 




i “nox j tlfpy, u, y J V^ — 2 m Y ~ X Tfiu 

m R 2 )n v h + 6 m R 2 f R<j) ( 0 , 0 , </> (0) ) = 

^ ^ 2 , , fM 0 , 0 , 0 (°)) 

(□r? + rriai )<p H-—- 


2 m R 2 XT 


The field equations (2.10) generalize those of Ref. [39], where no scalar field component were 
considered. Note also that these equations are the weak-field limit of the model discussed in 
Ref. [40], 

To solve Eq. (2.10) we introduce the auxiliary field 7 ^ such that 


(Dr; T m Y ) D ^ri'Y/iv 


(I— —177 T my )(T r/hf^u 

'm R 2 -m Y 2 , „ (™Y 2 , m R 2 + 2 ?ny 2 rn ^ 

3m R 2 ^ + M 2 + 6 m R 2 n J 

+2 my 2 f R<t> ( 0, 0, </> (0) ) (<9 2 v - rivvU^ip , 


D v h 

( 2 . 11 ) 


leading to 
h 


-/J,U — 7/71' T 

— 2 7715" 


m R 2 -m Y 2 2 /my 2 m fl 2 + 2 my 2 \1 2 , 

L 3 m * 2 + 2 + 6 m * 2 a v)\(Pn + ™ Y ) h 

2 /R 0 ( 0 , 0 , 0 (°)) (^ - ^□^(□„ + my 2 )- 1 n- : V . 


( 2 . 12 ) 


Since the trace h i 


, Wf) 

h = -£ 


is 

^(□^ + my 2 )(D ?7 + m R 2 ) _1 7-6mjj 2 / R 0(O,O,</) (o) )(n, ? + m R 2 ) _ V , (2.13) 

TYl y 


Equation (2.12) can be written 


as 


9 

7 m R 

" J fiv — Q - 


m 


Y L 


- 2 m A 2 /^( 0 , 0 ,/)) 


l 2 -m Y 2 2 , /my 2 mfi 2 + 2 m r 2 \1 2 x-i 

3^—^ + ^(~a~+ 6mj;2 + ^ 

n ^(°)y n-Un ™ ~ 2 "\—i,a 2 1 /'n 1 ™ 2^— 1 70 


Using Eqs. (2.12), (2.14), Eq. (2.10) reads 

(□»? + my 2 )D,7 F = -2 my 2 X T^ , 

(□, + my 2 )D, 7 = —2 m Y 2 X T , ( 2 - 15 ) 

(□r, + m^)ip - 3m R 2 f Rcj) (0, 0, 0 (o) ) (□,, + m, R )- l n v Lp = -m R 2 f R(j> { 0, 0, </> (0) ) (□,, + m* 2 ) _1 ;h T . 

Hence, Eqs. (2.10b) and (2.10c) have been decoupled. Let us rewrite Eq. (2.15c) as 

(□ 7 / + m+)(D v + m 2 _)ip = -m R 2 f R4 ,( 0, 0, (/> (0) ) X T , 
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(2.16) 























where 

m± = rn R 2 w± (£, rj) , 

with 

1 - £ + 'if ± \/(l - £ + r/ 2 ) 2 - 4?y 2 

™±(£,?7) = -g- 

£ = 3 /^ 0 ( O , 0 , 0 (o >) 2 , 
ms 

V = —» 
m R 

and the constraint £ < 1 . 

We then introduce the auxiliary fields T, T, H defined through 

(□»? + m R 2 )r = ~m R 7 , 

(□t? + m R )U v ^ = -2 m R 2 v, 

(I \3 v + m, R 2 )Z = -2 m R 2 ip, 


(2.17) 


(2.18) 


(2.19) 


so that the solution (2.14) reads 


hf!!; — + 2 


m 


Y L 


m-i? 2 - my 2 2 / my' 

Oui/ + I x I- 


3 m R 2 


m R 2 + 2 my 2 
Gm ^ 2 




+/^(O,O, 0 ( o) ) 


dlv * + ^Vuv 


( 2 . 20 ) 


In the limit m R —> oo, my —>• oo and for vanishing f R( /,(0, 0, (/>( 0 )) we recover the standard 
results of General Relativity (GR) since Eq. (2.20) reduces to h /M/ = 7 ^ 

To solve these equations we need the Green’s functions (see Appendix A); we hence 
consider the distributions C/KG.m, < 7 gr which satisfy the equations 

(□77 + m 2 )£ K G,m(ab x') = 5 a {x-x') , 

^ v Qgk{x,x') = 5 a (x-x') , (2.21) 

where l/KG,m and Ggr are the Green functions of a Klein-Gordon field with mass m (KG,m) 
and the standard massless models of standard General Relativity (GR), respectively. 

Due to causality, we are only interested in the retarded Green’s functions, hence 

T xx t ) / t 2 x , \ 

\ 2 / 

( 2 . 22 ) 

where 

T xx' = {x- x') 2 = (t- H f - |x - x '| 2 . (2.23) 

The terms with the Dirac distribution describe the dynamics on the light cone, i.e. t — t! = 
|x — x'|, while the ones with the Bessel function of the first kind J\fx) describe the dynamics 


Q&, m(x,x') = 

e r c?k(x,x') = 


Q^t — 


47r 

Oft— f 


S(t — if — |x — x'|) 


X — X' 


5{t — t' — |x — x'| 


47r 


X — X' 
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interior to the light cone, i.e. t — t! > |x — x'|. We can now build the Green’s functions for 
the auxilarly fields 7 ^, T, S and p as particular combinations of the Gkg m an d ^GR : 

a; et (x,x') = o( t -Ae(^f ) J ^ mYTxx ' ) , 

7 V / V 2 / 47 Tm Y Trr' 


e(t-n 0 


_„2 


XX' 

2 


G T r et (x,x') = 


Gt(x,x') = 


47r (m+ - mi) t X!E / 


0 (t - t ') 0 ( 


47 T (m 2 R - niy) t xx i 




47T7-, 


1 xx' 

m+ Ji(m+ t xx i) - m- Ji(m- t xx <) 
Ji (my Txx f ) _ Ji(m R x xx ’) 

my m R 

kii 1 (yn+ x X x ’) 


+ 


Ji(m r xx t) 


my (mi — ? 77 ,^_)(m^— m+) m_(my — mi)(m^,— mi) 


+ 


c7i (m/j ryy) 


mfi(mi — m\)(m_ — mi) 


£§*(*>®0 = - 


+ 


R) 

®{t - t'J 

m R Ji(m R 

7~xx', 


m + Ji(m + 

7XX ') 


+ 


T xx f) 


(mi — mi) (mi — mi) (mi — mi) (mi — m_) 


(m\ — m‘ R )(m 2 _ — m 2 R ) 

Using Eq. (2.24) we then derive the particular solution for the held 7 ^: 


(2.24) 


7^(i x ) = - 


my fh 
2 7T 


/ 


d 3 x' 

t— |x—X 1 


^ , J7i (myv/(t — t') 2 — |x — x'l 2 ) . A 

dt' n i_ : ; 1 =^T u Jt’, x') , 2.25 

■y/(f — t') 2 — |x — x/|2 " 1 ’ ^ ' 


where i — |x — x'| is the retarded time. Introducing the variable r = my W(f — t ') 2 — |x — x'| 2 , 
Eq. (2.25) takes the form 


7 /w(i x ) = - rr =!f X f dV I dr 


2t r 


/■ c 

C ' / 
do 


^i(t) 


r 2 + my |x — x '| 2 


x t - 


r 2 + m,y |x — x' 
my 


/12 


X 


in the limit my ->oowe obtain the standard results of General Relativity, namely 


7 y^(i 


x 


/ d :i x' 
2 7T ./ 


, - |x - x'|, x') 


X — X' 


(2.26) 


(2.27) 


3 Gravitational waves emitted by a quadrupole source 

Let us assume that the sources are localized in a limited portion of space within the neighbor¬ 
hood of the origin of the coordinates, namely the sources have a maximal spatial extension 
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I X maxI (where T^ v 7 ^ 0 if |x| < |x( Tiax |). If we consider the far zone limit, or radiation zone 
limit, i.e. |x| >> A >> |x' max | where A is the gravitational wavelength of the waves emitted, 
we can consider the solution (2.25) at a great distance from the source. In this limit (i.e. 
radiation zone), we set |x — x'| ~ |x| and the solution can be approximated by plane waves 
having nonzero only the spatial components, i.e 'ftt = 7 ti = 0 and 7 ^ 7 ^ 0. Note that in 
modified gravity models one has in general six different polarization states [39]. The spatial 
components of can be written as 


7 ij(t, x) » 



(3.1) 


where in general 


Tm. — 


\J t 2 + m 2 |x | 2 


m 


(3.2) 


The spatial components of the energy-momentum tensor Tij are related to the quadrupole 
moment 


Qij(t) = 3 J d 3 x T a (t, x')x'a= 3 J d 3 x p(t, x) x\ , 
through the relation 


1 d 2 


^x'Tj^t, x') = = 


Qij(t ) 


6 dt 2 ^' JS ~' 6 

Equation (3.1) can be casted in the form 

7ij( x ,t) = -2myT™ Y (lx\,t) , 


where 


X 


T™(t,\x\) = — I dr 


Ji(r) 


\Jt 2 + 


m 2 Ixp 


Qij Tm) 


Considering the trace of Eq. (2.25) in the radiation zone limit 

7 (t, x ) ~ ~ X [ dr = T ^ f d 3 x' T(t — r mY , x) , 

n \Jr 2 + rriy |x | 2 

we get 


J d 3 xT(t , x') = ff v J d 3 x T^ u (t, x',) = M 0 + 


Q(*) 


6 


(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 


where Mq is the mass of the source and Q(t ) = 'if' J Q VJ (t) is the trace of the quadrupole 
moment (3.3). Hence, Eq. (3.7) becomes 

,(*,<) = - ”*** r dr 


2 7T 
my T Z 100 

’ 12 7T J 0 


T 2 + |x | 2 


dr , ^ ^ ^ = Q(t - r my ) . 


r 2 + m,y | x | 2 


(3.9) 
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The first term on the r.h.s. of the equation above reads 

Jl(r) _ 1 


dr 


sj t 2 + m 2 |x | 2 


m x l 


1 — e 


(3.10) 


hence, it does not depend explicitly on time, so that it does not contribute to energy loss of 
the system. 

Using Eq. (3.6) we obtain 

Mr) 


t| x|) = ts, |x|) = A- 


24vr 


dr 


\J t 2 + m 2 |x | 2 


Q(t An) 5 


hence the solution for the trace 7 (see Eq. (3.9)) takes the form 

7(*> M) = —2 my T my (t, |x|) . 
The solutions for all fields then read 
7*i(t, l x l) = - 2 myT™ r (f, |x|) , 


(3.11) 

(3.12) 


7 (t, |x|) = -2m.yT“ y (f, |x| 


¥>(*, l x l) = - 


K-m^) 


?n + T m +(f, |x|) — m_T m -(t, |x| 


T(t, |x|) 


2m H 2 m r 2 T m v(t.|x|) _ T m H(t.[x|) 
(m' 2 R -m'^) m Y m R 


|x|) = 2m.ij 4 /^(0,0,^>(°)) 


T m + (Tjxj) 


m+ ( m^ — rri 2 , ) ( m 2 R — ? 


5) 


+ 


T m -(t,|x|) 

m_ (m5_— m^)( — m?_) 


(3.13) 


|x|) 


m#( rri 2 R )(rn 2 — m^) 


i(i, l x l) = -2m//^(0,0,()i (o) ) 


rn+T m +(t,|x|) 


m-T m - (i, |x|) 

“T 7772-™2 \/„,2-77T7 “T 


m fl T m fi(t.|x|) 


(m_ — m+)(m R — m+) (mq_— m_)( m R — m_) (m+ — m R )(m_ — m R ) 


After laborious mathematical calculations (see Appendix B for details) we can rewrite the 
spatial components of the perturbation h /w given in Eq. (2.20), as (from Eq. (B.5)) 


hij(\x.\,t) = -2my Y™j Y + r]ij\ m Y gy T my (|x|,t) - m R g R T mR (\x.\,t) 


E mS9S ^ ,r,) T ms (|x|,t) + gR 


2 

+ 3 


s=± 

^rd x M) 


B m y(|x|,t) B m «(|x|,t) 


my 


+ E 9s(t,ri) 
s=± 


my 

D% a i\*\,t) 


m R 


m s 


(3.14) 


This is the main result of our analysis which we will use in the following to constrain the free 
parameters of extended gravity models found in the literature. 
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4 Energy loss 


The rate of energy loss from a binary system, in the far-held limit, reads 


d£_ 

dt 


2 vr|x| 2 
5 A 


hijh ij . 


Let us make the choice 


Qij {t) = Qij cos(w(y) t + ^(ij)) + Qij , 
Q(t)= Qcos(yt + i?)+Q° , 


(4.1) 


(4.2) 


where Qij and its trace Q are the quadrupole oscillation amplitudes, Qd and its trace Q° 
are constant terms, topj\ and ’duj) are the frequencies of oscillations and the phases of the ij 
components respectively, while v and & are the frequency and phase of the trace, respectively. 
All these quantities are considered to be time independent. 

From Eqs. (3.14) and (4.1), the energy loss reads 


d£ 

dt 


m 2 Y ^ {ij) QijQ ij H 2 a 

720 vr 

i/ 6 Q 2 Ixl 2 X 




2880 vr 


^2 Cspmsmp F{\x\-,ms\m P \v) , 


{S,P} 


(4.3) 


where we have averaged over time neglecting higher order terms. Note that with the notation 
Xqs P} we ex t en d the sum over all possible values of {S', P} = {Y, R, +, —All quantities 
in Eq. (4.3) are defined in Appendix B. 

The model under consideration carries by itself a natural frequency scale linked to 
the mass scales mp, my, m+ and ?n_. The E(|x|;mi;m 2 ;o;) functions are highly oscillatory, 
with different behavior for uj > uj c and to < to c , while for to = to c are highly resonant [41], 
The to > to c case is excluded folowing a simple heuristic argument [42] we highlight below. 
A system with to > to c cannot decrease its orbital frequency across the lower boundary to c . 
Since one expects all astrophysical systems to have formed from the coalescence of relatively 
cold, slowly moving systems, it is reasonable to suppose that at some time in the past, all 
binary systems had to < to c . Hence, we will only analyze frequencies lower than to c . 

For to < to c the last function of Eq. (B.8), can be approximated by [41] 


F(|x|; mi; m 2 ; to) 


+ A yjm\ m- 2 |x|; 


mi m 2 x 


(4.4) 


where C ~ 0.175 is approximately constant except for to —> to c and A(x, y ) = C 2li- y 2 • 
From Eqs. Eq. (4.3), (4.4) we obtain a contribution from General Relativity (GR) and one 
from Extended Gravity Models (EGM): 


= —<^GR — £egm 
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d£ 

dt 


(4.5) 










where 


£gr = 


x 

720 7T 


,6 . 


uy^QijQv - 


u 6 Q 2 


(4.6) 


£egm = 


X 


720 7r W (ij)^ 


QnQ} j A 


my X 


U (ij) 


+ 


/ 6 Q 2 A” V Y,R,+,- 
2880 tt A^SP 


CspA ^/ms mp|x|; 


with 

(yy = 9Y (3 gy - 4) , Crp = 3 g R , Css = y , Cyr = ( 2 ~ 3g Y )g R , 

Cys = g (2-3 g Y )gs , Crs = gRgs , C± = ^g+g- , (4.7) 

for any values in {5, P} = {Y, R, +, —}. 

The correction to General Relativity, namely the term <?egMj ^ as ten characteristic 
frequencies (Table 1). We note however that oj+ + = uj r+ and = w^_, and in general 
u/j__ > u/j_ + > co ( L_ and Wy + > Wy_. Since the binary systems cannot have more frequencies 
than those predicted from the theory, it follows that a scalar-tensor fourth order gravity 
model has at most only five characteristic frequencies, i.e. Wyy, u c rr, uj yr , ui c __ and 
Note also if the trace of the quadrupole Eq. (4.2) does not depend on time then the correction 
(second line of Eq. (4.6)) depends only on Wyy. Therefore, theories constructed without the 
invariant R flu R^ u will not give different values than those of General Relativity. 


C^ TYly 

u c R _ = CrriRy/wZ 

Urr =Cm R 

w+_ = C m Ry /w + w- = C y/rriR 

w++ = Cm R yw^ 

Uy R = c y/my rn R 

UJ C __ = C VlRy/W- 

Uy™ c = C yjmy m Ry Jw^ 

to c R+ = C m R yw^ 

Wy = C yjmy mRy/W- 


Table 1. Ten characteristic frequencies for a scalar-tensor fourth order gravity. 


The correction for the energy loss given by an EGM model takes the form 


£egm = 


X 

720 7T 


Z2 6 Q 2 


"WjGyG* + 


Cyy 


A my |x|; 


U {ij) 


u>- 


YY 


+ 


v b Q 2 X 

2880 7T 


Crr A rtiR |x| ; —— + (c — + Cr-) A I 


m_ x ; 


u 


RR . 


5 C 

or 


+ Cyp A y/myrriR |x|; 


0J- 


YR, 


+ Cy- A y^my m_|x|; 


CJ 


Y— 


(4.8) 


As an example let us consider a pair of masses mi and m 2 in an elliptic binary system. For 
such a system, orbiting in the (x, y) plane, the nonzero components of the quadrupole (4.2) 
are 


On 

Qxx = 2 ^ ^ 


Qyy ~ 


3 l2 „ 3 3 n 9 

— 2 ® ’ Qxy — — /i a 0 , Q — — — fi a e , 


(4.9) 


with 


^xx — ^yy 


C ^xy — ^ — 2 ^ 


(4.10) 
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where /i = is the reduced mass, a and b are the major and minor semiaxis, e is the 

eccentricity and Cl is the orbital frequency. The energy loss then reads 


C _ O 6 f! 2 X 
t-GR — 20 ?r 


4(o 2 + 6 2 ) 2 —a 4 e 2 


c _ n 6 /j, 2 x 
‘--EGM — 20 it 


(4 (a 2 + b 2 ) 2 +a A e 2 (yy) A^my |x|; 


(4.11) 


+ a 4 e 2 (rr A \m R |x|; JA- + ((— + (r-) A m- |x|; £P- 


+ Cyr AI y/m Y m R |x|; ^ 1 +Cy-A[ ^Jm Y m_|x|; 


20 


These results will be used in the next section in order to constrain the three mass character¬ 
ising the extended gravity model under consideration. 

5 Observational constraints 

One has to test the observational compatibility of an extended gravity model. Hence we will 
study the variation of the orbital period V for binary systems due to emission of gravitational 
waves. The relation between the time variation of period and the energy loss is 


dV V 3 dS 


dt 47t/psr dt 1 


(5.1) 


where /psR is the pulsar’s moment of inertia, normally assumed to be 10 45 g cm 2 [43]. For an 
elliptic binary system in the weak-field limit, one gets 


where 


V 


4vt/psr 


(£gr + £egm) = T’gr + T*egm , 


(5.2) 


-rS _ Q 6 11 2 V 3 X 
' GR - — 80 7r 2 /psr 


4(a 2 + 6 2 ) 2 —a 4 e 2 


T*egm = — 


n 6 n 2 p 3 v 

80 tt 2 Ipsn 


4(a 2 + 6 2 ) 2 +a 4 e 2 ( YY 


A ( m Y |x|; £P- 


+ a 4 e 2 


(rr A ^RIr |x|; + (C—+ Cfl-) A^m_ |x|; 


(5.3) 


T Cy r A ^ 77i y m fl |x|; + ( Y _k^/m^ 


m_ x ; 


2 U 
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PSR J0348+0432 is a neutron star in a binary system with a white dwarf, with estimated 
massed (2.01 ± 0.04) M 0 and (0.172 ± 0.003) M 0 , respectively [43]. This binary system 
has an almost circular orbit, a semi major axis a = 8.3 x 10 8 m and a short orbital period 
V = 2.46 hours orbit. For these values, General Relativity leads to a significant orbital 
decay. In particular, the authors of Ref. [43] obtained the constraint 


Vobs/VGR = 1.05 ±0.18. 


(5.4) 


Using this result and Eq. (5.2), we get 


-0.13 < ' E g QM - < 0.23 

“ ^GR “ 


-0.13 < A 


s A my x 


20 


< 0.23 


hence 


my > 5 x 10 11 rn 1 . 


(5.5) 


Thus, for a binary system with a negligible circular orbit (e <C 1), one can always constrain 
the parameter my. To constrain the other parameters, one has to consider elliptic systems, 
i.e. the eccentricity must not be negligible. For example, for the elliptic binary system 
PSR B1913±16 [44, 45] where the experimental eccentricity is e = 0.6, the semi major axis 
a = 1.95 X 10 9 m, the orbital period V = 7.7 hours orbit, and Pobs/PcR = 0.997 ± 0.002, 
one infers 


-0.005 < 


4 (a 2 +& 2 ) +a 4 e 2 (yy 
4(a 2 +fe 2 ) — a 4 e 2 


A 


my x ; -c 


20 


+ 


4(a 2 +ft 2 ) —a 4 e 2 


Cl RR A 


x 


20 


± (C-+ Cr-) A 




+ Cyr A 


^y/my m R |x|; 



(5.6) 


+ Cy- A 


^y/my m— 



< - 0.001 . 


We have thus obtained a relation between the characteristic masses and frequencies for a 
general extended gravity model. In what follows we will examine some particular extended 
gravity models studied in the literature. 


6 Scalar-tensor fourth order gravity models 

Let us consider case A of Table 2; the only interesting quantity (see Eq. (2.8)) is rriR. For 
the system PSR B1913±16, Eq. (5.6) implies 

a 4 e 2 aIuir |x|; 

- 0.005 < -A- - L < -0.001 ^m R > 3 x lO'V" 1 . (6.1) 

- 4 ( a 2 + 5 2 ) 2 _ a 4 e 2 - 

In general, one can consider the polynomial expression 
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EGM 


Mass definition 


f(R, RapirP) 


™ 2 _ 1 
mR ~ 3/ m (0) 


my 2 —» oo, m^ 2 = 0 

£ = 0, Tj = 0 


m+ = m R , rri- = 0 


2 __1_ 

~ 3/ M (0,0)+2/y(0,0) 


my2 = TvTO’ m< ^ 2 = 0 


£ = 0, 1] = 0 


m + = m.R, m_ = 0 




/(-R, 0) + W(0)</>; Q </> ;C 


my 2 ->• oo, ?B 0 2 = -/^(O, ^»(°)) 

£ = 3/^(0, </>(°)) 2 , r, = % 


/ l-£+?r ! ± v '(l-f+r? 2 ) 2 -4 f? 2 

Wl± ~ V 2 - m R 


m R 2 ->• oo , my 2 -> oo, m/ = - 


D c 0 .R + ci R <£ + /(<£) + w(</>)</> ;Q </’ ;a £ = 3 ci 2 , i/ t 0 


m+ —> oo, m_ 


7 1—3ci 2 


^ _2 _ 1 

mR ~ ~ 3/fl fl (0, 0, ^(°))+2/ r (0, 0, ^(°)) 


my2 = / y (o,o, 0 (°)) ’ m <A 2 = -/^(°» °> ^ (0) ) 




£ = 3/ R0 (O,O>(°)) 2 , 17 = ^ 


/ i-C+^i^/fi-C+j? 2 ) 2 -^ 2 

m± — V - 2 - rn.R 


Table 2. Here /,r( 0, 0, </>^) = 1 and = 1/2 and for the case D we set also cq +Ci<// 0 ^ = 1 




( 6 . 2 ) 


N 

f(R) = R + aR 2 + UnR n - 

n = 3 

Note however that the characteristic scale niR is only generated by the i? 2 -term. An in¬ 
teresting model of f(R )-theories is that of Starobinsky f(R) = R — R 2 /Rq [46], for which 
m 2 j = Rq/6, hence using Eq. (6.1) we get Rq > 5.4 x 10 _19 m~ 2 . 

To generalize the previously result we must include the curvature invariant R pv R pu . 
For case B of Table 2 we consider the general class of f(R, R a pR a P )-theories and their 
characteristic scales rriR and my. Using Eqs. (5.5) and (5.6) we obtain 

my > 5 x 10 _11 m _1 , rriR > 1.15 x 10 _9 m _1 . (6.3) 

This class of theories includes the case of a Weyl square type model, i.e. C' Aa , p(7 C' /2; ' pfT = 
2R pi/ R p,u — | R 2 , where there is only one characteristic scale rriR —> oo. 

The same argumentation is also valid for the scalar-tensor case of theory, for which in 
the Newtonian limit (see case D in Table 2) the more general expression (2.6) becomes 

• (6-4) 

Thus, for the most general Scalar-Tensor (ST) theory in the Newtonian limit, one can consider 
the model 2 


/st (R, 4>) = Co R + c\ Rcj) — \m ct> 2 (0 - 0 (o) ) 2 + 


(6.5) 


Since for this case thr —> oo, my —> oo, £ = 3ci 2 , rj —> 0, m + —> oo and m_ 
we obtain from Eq. (5.6) 


^/l—3ci 2 ’ 


- 0.005 < 


a 4 e 2 


4(a 2 + 6 2 ) — a 4 e 2 


2ci 2 (3 - 6ci 2 ) 
(1-3 ci 2 ) 2 


nip |x| 2 Q, 

y/\ - 3 Ci 2 ’ 


< - 0.001 . 


( 6 . 6 ) 


As a special case of a scalar-tensor fourth order gravity model (case E) we consider 
NonCommutative Spectral Geometry (NCSG) [47, 48], for which at a cutoff scale (set as the 
Grand Unification scale) the purely gravitational part of the action coupled to the Higgs H 
reads [49] 


Sncsg = / 


R 


2 + c*o Cp Vpa C^P° + t 0 R*R* + - 


N ,H 

2 


R H 2 
12 


+ A 0 H J 


(6-7) 

where R*R* is the topological term that integrates to the Euler characteristic, hence nondy- 
namical. Since the square of the Weyl tensor can be expressed in terms of R 2 and R pv R pi/ 
as Cp U p (J C^ l ' prT = 2 R pu R pu — | R 2 , the NCSG action is a particular case of action (2.1). 


2 With the condition ao + an cf >^ = 1. 
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At this point, let us briefly discuss the behavior of the gravitational part of the action 
Eq. (6.7) above. We emphasise that Eq. (6.7) is the result of a perturbative expansion of the 
bosonic spectral action Tr(\'(Zl 2 /A 2 )), with % a cutoff action, D the Dirac operator and A a 
constant scale up to which the theory is valid; usually taken as the Grand Unified Theories 
scale [49]. The higher derivative terms that are quadratic in curvature can be written in the 
form [49, 50] 


d A x\/—g 


2rj 


C, 


/iispa 


C ^upa _ j± R 2 + 

3 r] 


E 


( 6 . 8 ) 


with E = R*R* the topological term which is the integrand in the Euler characteristic. 
Renormalisation Group Equations (RGE) determine the running of the coefficients g,uj,9 of 
the higher derivative terms in the action Eq. (6.8). The (weak) low energy constraints on the 
coefficients of the quadratic curvature terms R. IW R! W and R 2 impose that these coefficients 
should not exceed 10 74 [50], which is indeed the case for the action Eq. (6.8) as RGE analysis 
has shown [49]. 

Using Eqs. (5.5) and (5.6), we can constrain the parameter ao, which corresponds to a 
restriction on the particle physics at unification. We thus obtain 


«o < 10 20 m 


2 


(6.9) 


which is rather weak but can in principle be improved once further data of nearby pulsars 
are available. 

The parameter ao has been constrained in the past using either pulsar measurements [42] , 
Gravity Probe B or torsion balance experiments [51]. Here we have extended the original 
analysis of Ref. [42] for the case of pulsars with an elliptical orbit. Let us note that the 
strongest constraint on ao, namely ao < lCU 8 m 2 , was obtained [51] using torsion balance 
measurements. 


7 Conclusions 


In the context of extended gravity models and in particular within the class of scalar-tensor 
fourth order gravity, we have studied the energy loss of stellar binary systems. In general, 
the models we have considered depend on four parameters, as 


f(R,R a pR a P,<f>) + u;(#M> ; “ = 



2m fl 2 + my 2 r2 + R a/ 3R al3 
6mR 2 rriY 2 my 2 


(7.1) 





where we have set fn{ 0, 0, q i>(°)) = 1 and u (^(°)) = 1/2, Note that m^ 1 , rn y X , m^ 1 
determine the characteristic lenghts of propagation, £ is the coupling constant between the 
background geometry and the scalar field, while the background value <p 0 ' of the scalar field 
does not determine directly the dynamics. 

To constrain the free parameters we considered the nearly circular binary system PSR 
J0348+0432 and found my > 5 x 10 _11 m _1 . Considering the elliptic binary system PSR 
B1913+16 with eccentricity e = 0.6 we have constrained all three free parameters. Choosing 
a particular scalar-tensor fourth order gravity model (7.1) we obtained a lower value for m r 
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and m^, namely m R > 3 x 10 -9 m _1 and m > 1.3 x 10 -11 m _1 , respectively. One may be 
able to set stronger constraints by considering systems which are closer. It is worth noting 
that for circular binary systems there are no corrections in the case of a pure f(R) gravity 
with respect to General Relativity. 


A Green functions for a ScalarTensor Fourth Order Gravity 


The complete set of equations replacing the field equations (2.10) for tp and Eq. (2.19) 
for the auxiliarly fields T, T, S are 


(□»? + my 2 )D f? 7 M ^ = -2 m Y 2 X , 

(□, + my 2 )D ,7 = —2 my 2 XT , 

{U v + m+ 2 ){U ri + m- 2 )ip = -i V/R«/,(0.0,^ 0) )Tr, 

(A.l) 

(□,j + mij 2 )(D r? + ?ny 2 )D, ? r = 2 m R 2 m Y 2 X T , 

(n v + m + 2 )(n v + m. 2 )(n v + m R 2 )n v ^ = 2m R 4 f R<p (0 1 0 1 ^) XT , 

(D rt + m + 2 )(D v + m- 2 )(D 11 + m R 2 )E = 2 m R 4 f R(j) (0, 0, <^ (0) ) A T , 

where we have four characteristic lengths (m# -1 , my -1 , m+ _1 , m_ _1 ), which we assume 
all different and real. The first two ( m R ~ 4 , m Y _1 ) are generated by the geometry, while the 
last two (m_|_ _1 , m _~ x ) are lengths resulting from the interaction between geometry and the 
scalar field ip. 

The solutions of Eq. (A.l) can be expressed in terms of Green functions as 


7nv(x) = -2 m Y 2 X f d 4 x'G 1 (x,x')T llv {x') , 
j(x) = —2 m Y 2 X f d 4 x'(/j(x, x') T(x') , 

(p(x) = — m R 2 f R(j) (0, 0, (j)^) X f d 4 x'Qip(x, x') T(x') , 

(A.2) 

T(x) = 2m R 2 m Y 2 X f d 4 x'Gr(x, x') T(x') , 

T(x) = 2m fi %(0,0,<ji (o) ) X f d 4 x'Qq,(x, x') T(x') , 

E(x) = 2m R %(0,0,^°)) A f d 4 x'Q~(x, x') T(x') , 
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with the Green functions fixed by 


(□,, + my 2 )D I] 5 7 ( x,x') = 5 4 (x - x') , 

(□, ? + m \)(□,,+ m 2 _)Q^{x,x!) = 5 4 (x - x') , 

(\3 V + mR 2 )(n v + mY 2 )O n Qr(x,x') = 5 4 (x - x') , (A.3) 

(□, ? + m 2 + )(n v + ?n 2 )(□,, + m^ 2 )\J ri Qq,(x,x') = 5 4 (x - x') , 

(□,, + m+)(D r? + m 2 )(□,, + mR 2 )£s(x, x') = <5 4 (.t - x') , 


where 5 4 (x — x 1 ) is a four-dimensional Dirac distribution in flat space-time. To find the 
analytical dependence of the Green functions it can be shown that in Fourier space they are 
linear combination of only Gkg,ui and £?gr> which satisfy the second order equations 


(□r, + rn 2 )GKG,m(x,x') = S 4 (x-x ') , 
O v G G r{x,x') = S 4 (x-x > ) , 


(A.4) 


with solutions 


0 t-t' 

G r A = -\r 


e t-t' 


GgL(x, *0 = 4rr 


8(t-t'- |x-x'|) _ t xx , 

|x-x'| T xx , U l 2 


S(t—t' — |x—x'|) 


where t 2 x , = (x — x') 2 = (t — t ') 2 — |x — x '| 2 
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Hence the Green’s functions Ql ^, Q T J X . Q V ^ A and <5| et are 


G^M = A 


G%iM-G%h,m Y M 


— P) (+ — +'\ P\ ( A_\ '^l( rn Y T xx') 

- U V ) \ 2 ) **™YT XX , 


G T f(x,x') = v 1 V 

\ •> / mr — 




©^t-t'j ©^ 

47r (m^_— to^) t xx / 


m + Ji(m. + t xx i) - m_ Ji(m_ r xx /) 


gf et (x ,®') = , 1 , 

J- v ’ ' mA—m,Xr 




471- (m^-m^r) t xx , 


Ji( m YT xx ,) Ji(m R T xx ,) 


my 


m R 


fret/ A _ Q gr( X ! X ') _ £ KG,m + ( X ’ X ') _ ^KG,m^ ( x ’ x ') _ ^KG.mjj 

t^vji (T) *£ 7 mil m?, m‘ 2 ,(m^ — m 2 i)(m' R —m 2 ,) m 2 _ (— m'l) (mT — m 2 _) m 2 R (m 2 — ml)(ml- m R ) 


©I t-t' ] ©( 


JiKv) 


+ 




+ 


J\{m R r xx ,) 


m+(m_—mj_)(m|j—ml) m_( ml— m_)( rn 2 R — m_) m R ( m^_ — m R )( m_ — 


GA(x,x') = 


GF G ,mJ*y) 


G ret 

y KG,n 


(x,x') 


(ml — ml_)( rrij,— ml) " T " (ml — m 2 _){ m%— m 2 _) (ml— m 2 R ){ m 2 _ — m R ) 


%G,rajj' 


©( t-t' ) ©| ^ 


vi + Jl{m + T xx ,) r xx t ) m R Ji{ni R t xx ,) 

(m'f — m+)(m' R — m+) (m^ — m’j_)(m R — m 2 _) (m)( — rn R )( m 2 _ — m R ) 


B Mathematical aspects of the metric components /t^- 


The spatial components of the perturbation /t Mi , (2.20) can be expressed as 


1 


hij(t, |x|) = -2myT™ y (t |x|) + -- 2 <9f- + 77 ^ H T my (t, |x|) - 77 ^ n T m «(f, |x|) 


+ E 


s=± 


9s(Z,v) 

3 m s 


3 my - 


2 dfj - rjij m 2 s T ms (t, |x|) , 


(B.l) 


(A.5) 
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where 


9 m|jmy+3 (2 m^+my )□ 
(WR-m^) 

(2m|j+?ray)(n+m|j) 
mR ( m fl —m y) 


5± (£,d) 


_I_ 

[™TKfl)-»±Kfl)] [l—UI±(^,r?)] 


(B.2) 


and x) and T m (t, x) are defined in Eqs. (3.6), (3.11). The derivatives of T m (t, x) are 


^T m (t, |x|) 

3L/r m (t, |x|) 


|x|) 


□T m (t, |x|) 


X 

247T 


/•OO 

/ dr Ji (r) 

do 


m 2 ^(|x| 2 )Q(t - r m ) dfj,Q(t - r m ) 




m* x 


*y roc 

dTjlW 


m 
2 

v /*oo 

2 s/ dTjl(T) 

mQ'(t - T m ) 


2 [r 2 + m 2 |x| 2 ] 3 ^ 2 

3 m 4 9 At (|x| 2 )^(|x| 2 )Q(f - T m ) _ »r 2 ^(|xl 2 )Q(t - r m ) 

2 [t 2 + m 2 |x| 2 ] 3 ^ 2 
<9LQ(i - r m ) 1 


4 [r 2 + m 2 |x| 2 ] 5 ^ 2 

2 r ^(|x| 2 ) d v Q{t - Tm) + <9^(|x| 2 ) dfj,Q(t - Trr 


[r 2 + m 2 |x| 2 ] 3 ^ 2 

m 2 Q(t - r m ) 

[r 2 + ?n 2 |x| 2 ] 3/2 L[r 2 +m 2 |x| 2 ] 


+ 


a/t 2 + m 2 |x| 2 


3 ?7l 2 Xi Xj 


- Si 


+ 


2 m 2 x,- x 


* ^3 


T 

24vr 


[t 2 + m 2 |x| 2 ] [ [r 2 + m 2 |x| 2 ] 
. ...2 r 

dr Ji (r) 


-Si 


+ 


m 2 Xi Xj Q"(t — r m ) 

[r 2 + ?n 2 |x| 2 ] 


3/2 


3 771 2 Q(t- Tm ) T 2 + 3 m Q'{t - Tm ) T 2 


m 2 |X| 2 Q"(t - Tm ) 
[r 2 + m 2 |x| 2 ] 3 ^ 2 


[t 2 + m 2 |x| 2 ] 5 ^ 2 [r 2 + m 2 |x| 2 ] 

Q id T m) 


+ 


\J T 2 + m 2 |x| 2 


(B.3) 


Equations (B.3c) and (B.3d) can be approximated by considering only the terms scaling as 
l/|x|; the other terms scale as l/|x| n with n > 1. Thus, we have 


x2 ymu |„|\ „ jT f 00 ^ T-, f T ) d 2( 3(t Tm) 
v 1 ’ 1 |j 247r Jo ar ^l[ T J ^ T 2 +m2 | x |2 


24?r , ;g ^ D’Jdxl.t) 

T^+m^lxp 


□T 7 


■(«. |X|) « & r irJi (t) ^ g-(N, t) , 


(B .4) 
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and Eq. (B.l) reads 


| x |) = -2myT“ y + jmy g Y T my (t, |x|) - m R g R T nR (t, |x|) 


-Es=± mS9 f’ V) Tms (^ |x|)+5* 


B m Y(\xlt) B rn R{ |x|,t) 

m Y m R 


(B.5) 


+ 2 

' 3 


D™ Y (t, |*|) 


my 


+ J2s=± 9s(Z,ri) 


D™ s (t, |x|) 


where 


fly = 


3m|j 


2 2 
m R — my 


9R = 


m s 


2 rn 2 R + m-v 


y 


2 2 
171 R ~ m Y 


(B.6) 


The time derivatives of T™(t, |x|), B m (t , |x|), D^(t, |x|) needed to calculate the energy loss 
in Eq. (4.1) are 


|x|) = 


T m (t, |x|) = 


U (ij) Qij X 


24 TT 


v 3 QX 
24 7T 


sin ( w fo ’) 1 + %))/i ( m l x l; -£r - cos ( w fe) t + %))/i ( m l x l; 


"(O) 

m 


s 1 m Ixl; — 
l l ’ m 


S m (t, |x|) = - 


A?(*> l x l) = - 


24 7T 


sin(z/t + i?)/f ^ra|x|; - cos(z /1 + $)/f 


sin(z/t + $)/f^m|x|; - cos(ut + d)ff fm|x|; 


(B.7) 


z/ 5 m 2 aii ai 7 Q ^ 


24 7T 


sin 


(z^+ i?)/ 3 c ^m|x|; ^ — cos(z't + r?)/| ^ 


m x ; — 

1 I’m 


with the definitions 


\ roo , J'i (r) cos(z Vt 2 +x 2 ) 

f n (x;z) = f Q dr , 

„/ \ roo , ^i( T ) sinWT 2 +i 2 ) 

fn^ Z )=fo ^ ’ {v ^y , 


(B.8) 


F(|x|;mi;m 2 ;w) = /i (™i|x|; ^ /f fm 2 |x|; ^ + /f (mrM; ^ /f ^m 2 |x|; ^ . 
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